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Abstract We present the results on precision orbit deter- 
mination from the radio science investigation of the Lunar 
Reconnaissance Orbiter (LRO) spacecraft. We describe the 
data, modeling and methods used to achieve position knowl- 
edge several times better than the required 50-100 m (in total 
position), over the period from 13 July 2009 to 31 January 
2011. In addition to the near-continuous radiometric track- 
ing data, we include altimetric data from the Lunar Orbiter 
Laser Altimeter (LOLA) in the form of crossover measure- 
ments, and show that they strongly improve the accuracy of 
the orbit reconstruction (total position overlap differences 
decrease from ~70m to ~23m). To refine the spacecraft 
trajectory further, we develop a lunar gravity field by com- 
bining the newly acquired LRO data with the historical data. 
The reprocessing of the spacecraft trajectory with that model 
shows significantly increased accuracy (~20 m with only the 
radiometric data, and ~ 14 m with the addition of the altimet- 
ric crossovers). LOLA topographic maps and calibration data 
from the Lunar Reconnaissance Orbiter Camera were used 
to supplement the results of the overlap analysis and demon- 
strate the trajectory accuracy. 


Electronic supplementary material The online version of this 
article (doi: 10. 1007/s00 190-01 1-0509-4) contains supplementary 
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1 Introduction 

1.1 Objectives 

The Lunar Reconnaissance Orbiter (LRO) mission (Chin 
et al. 2007) was designed to provide a comprehensive (global) 
and detailed (high-resolution) survey of the Moon, in antic- 
ipation of the needs of future orbital and landed exploration 
missions. The payload consists of seven instruments (Von- 
drak et al. 20 1 0), some of which have very high spatial resolu- 
tion: better than 50 cm for the Lunar Reconnaissance Orbiter 
Camera (LROC) system (Robinson et al. 2010), ~5-m foot- 
prints and 10-cm vertical precision for the Lunar Orbiter 
Laser Altimeter (LOLA, Smith et al. 2010a); 320 x 160 m 
field-of-view for the Diviner Lunar Radiometer Experiment 
(DLRE, Paige et al. 20 1 0); 30-m zoom resolution for the Min- 
iature Radio Frequency Technology Demonstration (Mini- 
RF, Nozette et al. 2010). The LRO datasets will ultimately be 
combined in a single reference frame based on the LOLA- 
defined geodetic grid, which requires accurate and consis- 
tent trajectories for geolocation. As such, the LRO mission 
defined a position knowledge requirement: 50-100 m total 
position and 1-m radially (Vondrak et al. 2010). 

Our goal is to report on the current status of the LRO 
precision orbit determination (POD). To determine precision 
orbits, the LRO trajectory from 13 July 2009 to 31 January 
2011 was reconstructed continuously, except for short gaps 
near spacecraft maneuvers. As described below, the quality 
of the orbits is assessed primarily through overlap analysis, 
i.e., the consistency of trajectory segments when processed 
in two consecutive time periods. 
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To obtain the most precise orbits for LRO, we have devel- 
oped a gravity model based on the historical tracking data 
(Mazarico et al. 2010) and further tuned with data from LRO. 
This intermediary model is not optimized for geophysical 
purposes and is intended to be used only to improve the LRO 
position knowledge. The upcoming Gravity Recovery and 
Interior Laboratory (GRAIL) mission (Zuber et al. 2011), 
scheduled for launch in September 2011, should provide a 
lunar gravity field with far greater accuracy than what is pos- 
sible to achieve from analysis of LRO tracking data. 

1 .2 Previous work 

In the recent era of lunar exploration, the Clementine space- 
craft was the first to provide the data to determine the Moon’s 
gravity field to relatively high spatial resolutions (/ = 70, 
Lemoine et al. 1997). The subsequent Lunar Prospector, 
with its lower orbit, in particular during the extended mis- 
sion (<50km), enabled features as small as ~40 km (/ = 
150— 165) to be resolved (Konopliv et al. 2001). The recent 
JAXA SELENE mission consisted of three spacecraft: the 
main orbiter Kaguya, in a 100 km-altitude circular orbit, 
and two sub-satellites which participated in a VLBI exper- 
iment. One of those sub-satellites also provided relay capa- 
bility for the radio signals, between ground stations and the 
main orbiter. This enabled radiometric tracking of the far- 
side for the first time and led to great improvements in the 
determination of the lunar farside gravity field (Namiki 2009; 
Matsumoto et al. 2010; Goossens et al. 2011). 

1.3 Overview 

After a brief summary of the LRO mission as it relates to 
the presented work, we describe the tracking data types used 
here: radiometric Doppler and range, and altimetric cross- 
overs from LOLA. The orbit determination process is then 
described in detail, followed by a discussion of the various 
steps in the current-best orbit reconstruction: radiometric- 
only orbits, the addition of the altimeter data, and the inver- 
sion of a new LRO-tuned gravity field (LLGM- 1 , LRO Lunar 
Gravity Model), which is used to converge the final trajectory. 
Finally, orbit overlap analyses are supplemented by results 
using independent data: LOLA topographic maps of the polar 
regions and LROC calibration information. 

2 Data 

2.1 Mission summary 

The LRO spacecraft was launched at 22:32 UTC on 18 June 
2009 from Cape Canaveral Launch Complex 4 1 , and entered 
lunar orbit on 23 June 2009. Five maneuvers designed to 
gradually circularize the initial eccentric orbit followed, and 


the spacecraft commissioning phase was initiated on 27 June 
2009. The commissioning orbit was a quasi-frozen ~30 x 
200 km polar orbit (i ~ 90°), with its periapsis near the lunar 
south pole. 

The instrument commissioning phase began on 3 July 
2009. LOLA was initially turned on for 2 days at that time 
and started continuous data collection on 1 3 July 2009. Since 
then, it has operated continuously except for the monthly sta- 
tion-keeping (SK) maneuvers, and a February 2010 safehold 
event. 

On 15 September 2009, the spacecraft transitioned to its 
nominal 2-h period mapping orbit, with lower eccentricity 
and an average altitude of 50 km. Because of gravitational 
perturbations, the orbit eccentricity and argument of periap- 
sis naturally evolve, and monthly station-keeping maneuvers 
are necessary to preserve the mapping orbit (Houghton et al. 
2007). The eccentricity is 0.0054 ± 0.0019, and never more 
than 0.010, with the spacecraft altitude generally between 35 
and 65 km (referenced to a 1,737.4 km sphere). After 1 year 
of the nominal Exploration mission, dedicated to the survey 
of the polar regions and of 50 sites selected for NASA’s Con- 
stellation human space flight program, the LRO mission was 
extended for 2 years by the NASA Science Mission Direc- 
torate and began its Science mission. 

In addition to the monthly SK maneuvers, smaller angular 
momentum desaturation events (fill ) generally occur twice 
a month: once just before an SK maneuver, and the other 
about 2 weeks later. When the Sun-orbit geometry is favor- 
able (large /) angle 1 ), enough power can be generated with 
the solar panel held fixed, which reduces the angular momen- 
tum accumulated in the reaction wheels and can allow a SH 
to be skipped (e.g., in July 2010). 

Other kinds of maneuvers are less frequent (yaw flip every 
6 months, at ft ~ 0° to keep the single solar panel oriented 
towards the Sun), or unique (phasing maneuver to observe 
the LCROSS impact in October 2009). Table SI lists the 
maneuvers of the nominal and science missions. 

The orbital maneuver plan leaves LRO in free flight for 
~ 14-day periods, which is desirable for orbit determination 
and gravity field estimation. It enables orbit reconstruction 
to be executed over extended time spans, thereby facilitat- 
ing the recovery of long-period perturbations due to the long 
wavelengths of the gravity field (i.e., low-degree spherical 
harmonic expansion coefficients, especially the zonal and 
resonant terms). It also improves the contribution of altimet- 
ric crossovers to the solution, by reducing the number of 
effectively independent spacecraft initial states that need to 
be estimated. 

The regularity of the spacecraft SK maneuvers leads to 
a natural partition of the mission into monthly periods, 


1 The beta angle (ft) is the viewing angle of the Sun from the orbit, 
i.e., the angle between the orbit plane and the spacecraft-Sun vector. 
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or “phases”, each lasting for approximately one lunation 
(28 days). In the time period presented here (13 July 2009 
to 31 January 2011), there are 21 phases in total. Commis- 
sioning, which lacks SK maneuvers, is divided into three 
phases of equivalent duration (CO_01 to CO_03). The nom- 
inal mission is split into 13 phases (NO_01 to NO_13), and 
the currently ongoing science mission consists of 5 phases 
(SM_01 to SM_05), as of February 2011. Table 1 gives the 
temporal extent of each phase. 

2.2 Data description 

2.2 . 1 Tracking data 

Radiometric tracking is the most widely used method of 
tracking interplanetary spacecraft, providing good measure- 
ment accuracy and high operational availability. As for the 
case of the Lunar Prospector (LP) spacecraft in 1998-1999, 
the LRO spacecraft is tracked continuously, i.e., anytime it 
is visible from one of the ground network stations. However, 
unlike LP which was tracked by the NASA Deep Space Net- 
work (DSN), LRO is tracked by a new LRO-dedicated NASA 
station in White Sands, New Mexico, and by the commercial 
Universal Space Network (USN). In Table S2, we list the sta- 
tion positions and their contributions to total tracking. All of 
those stations produce S-band Doppler and range radiomet- 
ric data, but only the White Sands station handles the science 
data downlink (telemetry), at Ka-band frequencies. The data 
are corrected for tropospheric delays from meteorological 
data collected at the stations. 

The tracking strategy and S-band tracking precision 
requirements (Imms -1 for the White Sands station and 
1.5-3 mm s -1 for the USN stations) were chosen because 
of the mission cost and the navigation positioning and pre- 
diction requirements (500-m total position), but were too 
large to enable the position reconstruction accuracy require- 
ments of the final LRO products (50-m total position). For 
that reason, a Laser Ranging capability was added to the 
spacecraft, providing one-way range measurements with 
a precision of 5-10 cm averaged over 5 s (Zuber et al. 
2010). However, the radiometric tracking performance is 
better than anticipated (—0.3 mm s _ 1 and —0.2 m for White 
Sands, ~0.4-0.8mm s” 1 and ~0.4 m for USN; Table S2), 
and greatly improves the achievable orbit accuracy. For this 
reason and because of the intricacies of the Laser Ranging 
data (in particular, the estimation of station and spacecraft 
clock and time bias parameters), early POD efforts and the 
results from the current work only include radiometric track- 
ing data. The use of Laser Ranging in LRO orbit reconstruc- 
tion is however a subject of active research. 

The radiometric data are complicated by the presence of 
biases which are typically not present in DSN tracking data. 
The Doppler measurements are biased, which is problem- 


atic because they usually provide the strongest control of the 
spacecraft orbit (Table S2). Early in the mission, the Doppler 
data collected by the White Sands station were time-biased by 
1 /76th of a second (around —13.16 ms); this was corrected 
on 19 October 2009. The Range data from all stations are 
also affected by timing biases. Although these can be rather 
well characterized for each station (—5 ms for White Sands, 
and —2 to —4 ms for USN stations), they are not constant per 
se, and still need to be systematically estimated. Figure S 1 
shows the range residuals from a sample arc (converged), 
with and without timing bias estimation. 

2.2.2 Altimetric data 

When available, accurate altimetric measurements can aug- 
ment the quality of the orbit reconstructions. At each 
intersection between two altimetry tracks (a “crossover” 
location), a constraint linking the height measurements made 
at the two epochs and the surface height change can be cre- 
ated. In the case of the Moon, with small-amplitude tides 
(—5-10 cm; Williams et al. 2008), no physical topographical 
change above the instrument noise level (10 cm) is expected, 
so the reconstructed altitudes of the altimetric bounce points 
should be identical. 

Rowlands et al. (1999) used altimetric crossovers for the 
orbit reconstruction of the Mars Global Surveyor (MGS) 
spacecraft and instrument pointing offset estimation of the 
Mars Orbiter Laser Altimeter (MOLA; Zuber et al. 1992). 
A limited amount of MOLA crossover data (—21, 000 over 
a 25-30 day period) was used in the orbit determination 
of MGS and Mars gravity field inversion (Lemoine et al. 
2001). A one-month simulation study of LRO (Rowlands 
et al. 2009) indicated that significant improvements could be 
obtained from LOLA’s multi-beam configuration. LRO is the 
first NASA lunar spacecraft to carry a laser altimeter with suf- 
ficient capabilities to allow the use of altimetric crossovers in 
the orbit determination process. Goossens et al. (201 1) noted 
improvements in certain situations with the inclusion of al- 
timetric crossovers on the SELENE spacecraft, even though 
the 1-Hz LALT instrument (Araki et al. 2009) has poorer 
spatial resolution than the 28-Hz LOLA instrument. 

The Lunar Orbiter Laser Altimeter is a 10-cm-precision 
28-Hz, five-beam laser altimeter, described in detail by 
Smith et al. (2009). At the typical 50-km altitude, the five 
5-m-diameter footprints form a 25-m-radius cross, slanted by 
26 degrees with respect to the direction of flight. With each 
laser shot, this pattern is repeated 56 m downtrack, and pro- 
gressively produces five profiles, separated in the cross-track 
direction by 10 m. Early results were presented by Smith et al. 
(2010a). With its high pulse repetition rate and multi-beam 
capability, the instrument collects up to 140 measurements 
per second. As of 31 January 2011, the LOLA dataset con- 
sists of 3.42 billion valid measurements. 
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Fig. 1 a Time series of the ft angle (in degrees) between 1 3 July 2009 
and 31 January 2011. For reference, the<$// maneuvers are indicated by 
a dashed line and the other maneuvers by a solid line, b Orbit-averaged 
LOLA return rates (i.e., percentage of ground returns) for spot 1, spot 3 
and the 5 spots combined. Returns for spots 2 and 5 are similar to spot 
1 ; spot 4 follows spot 3 closely. The early lower return rates are due to 
the eccentric commissioning orbit, when LRO was flying high over the 
northern hemisphere. Later, the repeated dips are correlated with high 
ft angle values, when the spacecraft is always near the terminator, and 
are the result of the ‘LOLA anomaly’ 

Compared to the MOLA data on Mars (Zuber et al. 1992; 
Smith et al. 2001), with its single-beam 160-m-diameter 
footprints spaced by ~300 m, the higher resolution of the 
LOLA measurements improves the coverage around the 
groundtrack intersections, enabling more numerous (up to 25 
effective crossovers at each groundtrack intersection because 
of the five beams) and better-quality constraints through more 
reliable interpolation. 

However, on the nightside, the data collection is ham- 
pered by a thermal alignment anomaly (Smith et al. 2010a). 
When the instrument is exposed to cold temperatures, the 
thermal blanket contracts pulling the laser beam expander 
out of alignment with the receiver telescope. The displace- 
ment is significant (~600 |irad) and the laser spots move out 
of the detectors’ fields of view (~ 400 |xrad in diameter). For- 
tuitously, the displacement magnitude and direction (along 
the Y-axis in the LOLA frame) combine to pull the transmit- 
ted laser beams 2 and 5 into the fields of view 3 and 4, so that 
altimetric data can still be collected the majority of the time. 
Consequently, channels 3 and 4 show much better average 
return rate than the others, as illustrated in Fig. 1. 

Every orbit, altimetric data cannot be collected over short 
periods (typically less than a minute) around the termina- 
tor crossings, i.e., at the poles. Although the data density is 
reduced, these important regions can still be well mapped 
thanks to the polar orbit convergence and changing illumi- 


nation conditions with season. In periods of high /I angle 
(|/3 1 ~ 90°), when the LRO orbit is nearly aligned with the 
terminator, the orbit-averaged lunar surface temperature is 
low, and the orbit-averaged return rate (percentage of good 
returns) drops dramatically (Fig. 1). However, the depen- 
dence on the ft angle is very steep, and the instrument shows 
substantial instrument performance degradation only a few 
weeks per year. 

The m-daily perturbations of the lunar geopotential (Kaula 
1966) produce periodic perturbations in inclination; the most 
significant terms are the m = 1 to m = 4, with periods of 
28, 14, ~10, and 7 days, respectively. Over the course of a 
month, the peak-to-peak variations are ~1.3°. These incli- 
nation changes, on an otherwise polar orbit, allow ground- 
tracks to intersect and altimetric crossovers to be computed. 
Between 13 July 2009 and 31 January 2011, we found 
5,007,658 crossover locations, potentially a slight overes- 
timate given our approach (valid data segments are picked 
liberally to avoid missing any crossover). When the data 
are actually prepared, only crossovers with sufficient valid 
LOLA data are used to ensure the quality of the interpolation. 

Figure 2a shows the distribution of those crossovers. It is 
very heterogeneous, with extreme polar dominance, due to 
the polar convergence of the LRO polar orbit: only 12, 5.7 
and 3% of the crossovers occur at latitudes less than 80° , 60° 
and 40°, respectively. The drape-like pattern is the result of 
the slow rotation of the Moon and the LRO high-inclination 
orbit, as noted by Rowlands et al. (2009), and the four longi- 
tude bands with fewer crossovers are due to the SK maneu- 
vers, which always occur at the same location, and lead to 
the repeat of those monthly patterns. Figure 2b is identical 
to Fig. 2a, but with a resolution of 0.25° pixel -1 instead of 
2° pixel -1 , which illustrates the fact that while saturating the 
polar regions, all the available LOLA crossovers so far only 
sparsely cover the equatorial regions. 

There are three main consequences of the LOLA anom- 
aly for the POD work presented here. First, the periods near 
ft = 90° contain fewer crossovers because of the lower aver- 
age return rate. This can affect the orbit reconstruction qual- 
ity of certain phases, in particular NO_03, NO_10, NO_l 1, 
and SM_05. In the future, this problem could potentially be 
mitigated by the use of multi-month crossover constraints 
(discussed in Sect. 3.3). 

Second, the general strength of the crossovers is reduced. 
Because of the slow rotation of the Moon and the polar orbit 
of LRO, the crossovers occur predominantly between ascend- 
ing and descending tracks. Given that we process crossover 
data in monthly batches, this means that in most intersecting 
track pairs one will have been acquired on the dayside, and 
the other on the nightside. Thus, instead of the anticipated 
geometry of five profiles intersecting five others (Rowlands 
et al. 2009; Mazarico et al. 2010), we usually have two tracks 
intersecting five and only ten single-beam crossovers (the 


■£) Springer 


E. Mazarico et al. 



Fig. 2 Maps showing the counts of LOLA groundtrack intersections in 
2x2° (a) and 0.25 x 0.25° (b) bins. The color scales focus on the lower 
counts to show the non-polar distribution. The pole values are strongly 
clipped, with actual count maxima of 17,929 (a) and 1,256 (b). The 
total count for (a) and (b) is 5,007,658. Frames (c) and (d) show simi- 
lar maps, but of the counts of the LOLA crossover locations that were 


actually used in the orbit detennination process (i.e., only the cross- 
overs occurring within each phase). The clipping is not as pronounced, 
although the maxima are 1,276 (c) and 200 (d). The total count for (c) 
and (d) is 172,293. Frames (b) and (d) show that, although the coverage 
is global, the possible sampling and actual sampling are still sparse 


measurement type used here, similar to the MOLA cross- 
overs used by Lemoine et al. 2001) are available, on average, 
at each crossover location. 

Third, for the same reason, much less “cross-track” infor- 
mation is available in one of the two intersecting tracks, 
which significantly interferes with the recoverability of the 
relative geometric adjustment of both tracks. Only the (sin- 
gle-beam) altimetric crossovers of the type used by Lemoine 
et al. (2001) with MGS/MOLA are employed, albeit more 
comprehensively than for MOLA and over a longer time 
period. 

3 Modeling and analysis 

3.1 Modeling 

Precision orbit determination of the LRO spacecraft is per- 
formed with the GEODYN system of programs developed 
at NASA Goddard Space Flight Center (Pavlis et al. 2006). 
The spacecraft trajectory is integrated using a priori force 
models, with a number of model parameters estimated iter- 
atively through batch least squares from the residuals of the 
observations with respect to the computed (modeled) values. 


GEODYN was used in previous POD and gravity inver- 
sion work for the Moon (Lemoine et al. 1 997 ; Mazarico et al. 
2010; Matsumoto et al. 2010; Goossens et al. 2011), Mars 
(Rowlands et al. 1999; Lemoine et al. 2001) and Mercury 
(Smith et al. 2010b). Most of the supporting data, and force 
and measurement models are similar to those described in 
previous studies. The lunar orientation and ephemeris are 
based on the JPL DE421 ephemeris (Williams et al. 2008). 
We used the degree 150 GLGM-3 gravity field (Mazarico 
et al. 2010) as a priori. That gravity solution is the result of 
the reprocessing at NASA GSFC of historical lunar tracking 
data (Lunar Orbiters, Apollo sub-satellites, Clementine and 
Lunar Prospector), in preparation for the LRO orbit deter- 
mination effort. The purpose was to create normal equations 
of the existing data for integration with the upcoming LRO 
data, as described in Sect. 3.3. The a priori value of the Love 
number ki is 0.027 (following Lemoine et al. 1997 who used 
the value determined from Lunar Laser Ranging by Williams 
et al. 1987). The trajectory integration is performed with a 
5-s time step. To compute the non-conservative accelerations, 
the LRO spacecraft is modeled as a 10-plate macro-model 
and is oriented in space according to attitude telemetry in 
the form of quaternions. The self-shadowed plate areas are 
computed as described in Mazarico et al. (2009). The lunar 
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surface albedo model is based on a low-degree spherical har- 
monic expansion (Floberghagen et al. 1999). Planetary ther- 
mal radiation was not modeled, but its expected effect is 
small (a few meters, Mazarico et al. 2010) over the typi- 
cal arc length. Following Mazarico et al. (2010), a constant 
along-track empirical acceleration is estimated to reduce the 
along-track orbit reconstruction error; the adjusted accelera- 
tions are small, around —2 ±6 x 10 -l() m s -2 . The center- 
of-mass position with respect to the spacecraft body frame is 
considered variable and is computed based on the fuel mass 
and the orientation of the solar panel and high-gain antenna 
(HGA). Because the two gimbals of the LRO FIGA orienta- 
tion system are separated by ~ 10 cm (which is the LOLA-LR 
instrument precision), the HGA phase center offset is mod- 
eled by two non-juxtaposed single-axis gimbals. The laser 
altimeter data are modeled from the LOLA receiver telescope 
position, and constant roll and pitch pointing biases around 
the post-calibration boresight vector are estimated. The data 
weights were set to values several times above the intrinsic 
noise and residual root mean square (RMS): 1 mm s -1 for 
the radiometric Doppler, 10 m for the radiometric range, and 
1 m for the altimetric crossovers. 
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Fig. 3 Doppler (a) and range (b) residuals are plotted against the cosine 
of the a angle (Earth viewing angle, i.e., elevation of the Earth from the 
orbit plane). Face-on or edge-on (as seen from the Earth) geometries 
plot to the left or right, respectively. The residuals obtained after con- 
vergence with two a priori gravity fields (GLGM-3 and LLGM-1) are 
shown 


3.2 Initial data analysis 
3.2.1 Radiometric-only orbits 

We process the data in several stages because the measure- 
ment types require different levels of processing effort: the 
radiometric data can be used readily and independently, but 
the altimetric crossovers, although they strengthen the solu- 
tion, cannot stand on their own. 

The radiometric data are divided into time periods called 
‘arcs’, to be integrated and converged individually. Follow- 
ing previous work on Lunar Prospector (Konopliv et al. 2001 ; 
Mazarico et al. 2010), we used short arcs, with an aver- 
age duration of 2.5 days. This length allows the arc to start 
and end with tracking passes from the higher-quality White 
Sands station. Although not necessary, this 8- to 12-h over- 
lap period proves an important and consistent monitor of 
the orbit reconstruction quality, as discussed below. These 
short arcs are adapted to high-quality orbits: long enough 
to include three whole White Sands tracking sessions and a 
number of intervening USN tracking passes, but short enough 
to avoid the detrimental build-up of modeling-related errors. 
In total, we constructed 272 short arcs to fully cover the 21 
analyzed phases. Figure S2 shows a typical arc distribution 
over 1 month. 

We process the data in monthly batches, well after they 
have been collected. Each phase is divided into 12 to 14 short 
arcs, which are individually cleaned and converged. In this 
first step, we only use the radiometric data. The initial state 
(position and velocity) of the spacecraft is obtained from the 


navigation team (Flight Dynamics Facility, at NASA GSFC) 
daily reconstructions. After a one-iteration run to manually 
edit out the anomalous Doppler and Range measurements, 
the arc is iteratively converged, by jointly estimating the 
values of typical arc parameters (initial state, measurement 
biases, radiation pressure scale factor and empirical constant 
along-track acceleration). 

The data reduction is generally of good quality, with 
Doppler and range RMS values between 0.4 and 1 .5 mm s -1 
and between 1 and 4 m, respectively (Fig. 3). The arc RMS 
values are strongly correlated with the Earth viewing geom- 
etry angle a~: lower for face-on geometries (cos a = 0) and 
higher for edge-on (cos a = 1, when the orbit crosses the 
deep far side). The better fits do not indicate better orbits: 
when face-on, the line-of-sight direction is aligned with the 
horizontal component of the gravitational residual acceler- 
ation and expected to be smaller than the radial component 
observed in edge-on orbits. 

We can assess the orbit precision (repeatability) by evalu- 
ating the RMS position difference of the trajectories com- 
puted in consecutive arc pairs during their overlapping 
period. This is done in the along-track (A), cross-track (C) 
and radial (R) directions (ACR), as well as in terms of total 
position (T). Over the whole period, there are 224 overlaps. 
(Some potential overlapping periods are voided by maneu- 
vers.) The overlap duration is 10.5 ± 2 h (see Figure S3 for 
the distribution histogram). 


2 The alpha angle (a) is the viewing angle of the Earth from the orbit, 
i.e., the angle between the orbit plane and the spacecraft-Earth vector. 
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Fig. 4 Phase-average overlap 
RMS values in the along-track 
(a), cross-track (b) and radial (c) 
directions and in total position 
(d) for various gravity fields and 
data usage. Blue indicates 
GLGM-3 radiometric-only 
orbits; cyan GLGM-3 with 
radiometric and altimetric data. 
Similarly, LLGM- 1 
radiometric-only orbits are 
shown in red, and the orbits with 
additional altimetric crossovers 
in magenta 
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The quality of the orbits reconstructed with the radiomet- 
ric data alone is close to the LRO requirements, although 
some periods can show larger discrepancies. In Fig. 4, we 
show the “phase”-averaged overlap RMS time series. Table 3 
gives an overall summary of the overlap differences; 47.6 m 
along-track, 49.9 m cross-track, 4.5 m radial, and 70.1m total 
position. The overlap RMS values calculated over each phase 
are presented in Table S3. In general, along-track and cross- 
track overlaps are smaller in edge-on and face-on geometries, 
respectively (Figure S4). 

3.2.2 Addition of altimetric crossovers 

Once all the arcs of a given phase have been converged with 
the radiometric data, they are ready to be combined, which 


is a pre-requisite for the use of altimetric crossovers. Indeed, 
given the slow rotation of the Moon and the polar inclination 
of LRO, altimetric crossovers mostly occur between tracks 
~14 days apart (or any multiple). 

One month of individual short arcs are processed simulta- 
neously by GEODYN, similar to the analysis of GPS data 
(Luthcke et al. 2003). Altimetric crossovers provide the 
“connection” between the separate data arcs. The ~28-day 
duration is chosen to allow any orbit to share an altimet- 
ric crossover with another orbit, which would not be the case 
with a 14-day duration: only the very first and very last orbits 
would have crossovers. Crossover computations inevitably 
require non-sequential accessing of orbital information. On 
each iteration, the numerical integration of orbits, and asso- 
ciated partial derivatives, is carried out for the whole arc in 
advance of crossover computations, and only certain quan- 
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Table 2 Comparison of the values obtained by various studies for key low-degree gravitational coefficients: GM (a), offset by —4.9028 x 
10 12 m 3 s -2 ; C20 (b), offset by +9.08938 x 10~ 5 ; C21 (c); S21 (d); and £20 (e) 


Reference 

GM - 4.9028 x 
10 12 (m 3 s“ 2 ) 

(C 20 + 9.08938 x 
10^ 5 ) x 10 9 

C 21 x 10 9 

5 2 i x 10 9 

£20 

Dickey et al. (1994) 

- 

- 

- 

- 

0.03020 ± 0.0012 

Williams et al. (2009) 

- 

- 

- 

- 

0.0199 ± 0.0025 

Williams et al. (2011) 

- 

- 

- 

- 

0.0229 ± 0.0020 

LP100K 

238,000 + 20,600 

20.009 + 5.319 

8.383 + 4.292 

7.707 ± 0.350 

n/a 

LP150Q 

1,076,100 + 8,100 

-7.295 ± 5.348 

-1.863 + 4.409 

-1.425 + 0.335 

0.0248 ± 0.0030 

Goosens and Matsumoto 2008 

- 

- 

- 

- 

0.0244 ± 0.0080 

SGM100H 

2,121,666+ 13,000 

-91.596+ 1.631 

11.270+ 1.064 

-15.894+ 1.443 

0.0240 + 0.0015 

SGM100I 

225,153 

-102.97+ 1.650 

7.861 ± 1.045 

-10.809+ 1.349 

0.0255 ± 0.0016 

SGM150 

795,109 

6.832 + 2.463 

28.776+ 1.534 

8.248+ 1.547 

n/a 

GLGM-3 

238,000 (not adjusted, 
=LP100K) 

3.654 + 5.344 

-29.620 + 7.719 

27.715 ±-4.357 

0.0235 ± 0.0020 

LLGM-1 

105,594 + 9,980 

5.102 + 1.916 

-52.865 + 1.330 

-13.345 + 1.570 

0.0270 ± 0.0004 


References are given in the text. The uncertainties quoted are the 1 -sigma values, except for £20 where we report ten times the spacecraft-derived 
formal errors following Konopliv et al. (2001) and Matsumoto et al. (2010) 


tities from the integration are stored in memory. This way, 
a balance is struck between memory requirements and the 
computational requirements associated with non-sequential 
accessing of orbital information. 

Because of the LOLA anomaly, we know that the laser 
altimeter alignment changes significantly between the day- 
side and the nightside, and that it varies rapidly near the 
terminator. However, little data are collected during those 
short transition times (Sect. 2.2.2), which makes it difficult 
to characterize the alignment other than by constant offsets 
from the nominally calibrated pointing. Thus, two indepen- 
dent sets of altimeter pointing biases (in roll and pitch) are 
estimated: one for the five daytime spots, and another for the 
two nighttime spots. Because of the large disparity in cross- 
over temporal density, the adjusted values are really reliable 
only when estimated once a month (i.e., the same adjustment 
applied to all the arcs). 

With the multi-satellite approach outlined above, we are 
restricted to a narrow subset of the ~5 million LOLA cross- 
over locations at hand (Sect. 2.2.2). Because we process data 
in monthly batches, only the crossovers occurring within each 
month (“intra-month crossovers”) are available. Multi-month 
arcs, designed to access the more numerous “inter-month 
crossovers”, are possible but computationally challenging 
and of limited interest to orbit reconstruction work due to 
the diminishing return. Figure 2c, d shows the spatial dis- 
tribution of the crossovers actually used. For each phase, 
we choose the crossovers which have enough LOLA data 
(for interpolation quality) and yield good adjustments, and 
then select the appropriate individual range measurements. 
Table 1 indicates the number of crossover locations available 
within each phase, with a strong /3 dependence (Sect. 2.2.2). 


Outofthose~172, 000 crossover locations, about ~ 159, 000 
are retained by GEODYN and contribute to the trajectory 
adjustment. In terms of actual measurements used, ~95 1 ,000 
(single-beam) crossovers exist at those ~ 1 72.000 locations, 
of which ~37,000 are rejected during convergence. Editing 
criteria are high RMS of polynomial fit to the topography, 
high-slope terrain, large altitude discrepancy between the 
two tracks and sensitivity with respect to the arc parameters 
(i.e., large state partial derivatives). The dominant factor for 
LRO is the latter, with some crossovers inclined to disturb 
the trajectory too much compared to what we expect from 
the majority. Compared to the MOL A analysis (Lemoine 
et al. 2001), the maximum slope was relaxed (from 0.25° 
to 1.0°), thanks to the shorter along-track distance between 
consecutive bounce points and the resulting more reliable 
interpolation. 

After orbit adjustment, the remaining height differences 
(crossover discrepancies) are small and comparable to the 
radial direction overlap (Tables 1, 2). The overall residual 
RMS is ~2.4 m. As noted above, some months, in particular 
at high f3 angles, contain fewer crossovers than nominally 
possible. In those cases, the crossover measurements do not 
significantly help the orbital solution. Figure 5 shows the 
distributions of both edited and retained crossover discrepan- 
cies, for a strong phase (NO_01 ) and a weak phase (NO_03). 
Only ~12% of the crossovers are rejected in the first case, 
but ~49% of an already smaller number are discarded for 
NO_03. Nevertheless, the RMS discrepancy does not exceed 
5 m (cf. Table 1). As shown in Table S3, the phase NO_03 
deteriorated in along-track and radial consistency, and the 
crossovers provided only marginal total position improve- 
ments. This is in contrast with the spectacular enhancement 
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Crossover discrepancy histogram 



Crossover discrepancy (m) 

Fig. 5 Histogram of the crossover discrepancy (interpolated altitude 
difference at point of intersection of two LOLA spot groundtracks) for 
two different multi-satellite arcs: phase NO_01 (2009-09-26 to 2009- 
10-23, in blue) and phase NOJ)3 (2009-11-20 to 2009-12-17, in red). 
In each case, the solid line shows the distribution for the (single-beam) 
crossovers retained in the solution and the dashed line that of the rejected 
crossovers. While the strong NO_0 1 month has only ~ 1 2% of its numer- 
ous crossover measurements excluded, the NO_03 month, with already 
scarcer crossovers, is further weakened by the rejection of nearly half 
the crossover data. Note the logarithmic scale for the horizontal axis 

observed over most orbital phases (Table S3): the overall 
overlap RMS values decrease from ~70 m in total position 
to ~23 m (Table 3), a 70% reduction. 

Obviously then, the use of altimetric crossovers in the 
orbit determination process can help meet the LRO position 
knowledge requirements, also showing that they can substan- 
tially mitigate the gravity-related modeling errors. 

3.3 Gravity field inversion 

The natural next step is to use the signals still present in 
the measurement residuals to improve the gravitational force 
modeling. We use GEODYN to create one set of normal 
equations per phase and per measurement type (i.e., Dopp- 
ler, range and altimetric crossover separately, which allows 


re-weighting of the various measurement types during the 
inversion). For practical reasons, we worked only on phases 
CO_01 to NO_13. 

Substantial changes were made to GEODYN to enable the 
computation of the normal equations for the crossover mea- 
surements. We mentioned above the need to pre-in tegrate the 
trajectory to use these data for orbit convergence (i.e., store 
the partial derivatives of the initial state and other arc param- 
eters), but this need is exacerbated by the large number of 
gravitational Stokes coefficients to be estimated (22,797 for 
a degree-150 gravity field). We use a multi-step approach, 
with large temporary files to store the necessary information. 

In addition, while the overlapping short arcs are important 
in the assessment of the orbit quality, it would be misguided 
to have some of the data (every other White Sands pass) 
counted twice in the normal equations. Instead of just trun- 
cating most of the short arcs by ~12h (to remove the last 
White Sands pass, for example), which would make them 
weaker, we created an alternate set of 91 longer ~5-day arcs 
(compared to 207 short arcs over the same time period). The 
quality of those arcs can of course not be assessed directly 
through overlap analysis, but orbit differences with the bet- 
ter-characterized short arcs are small. The average RMS orbit 
differences are 26.5, 13.5, 2.8 m in the ACR directions, and 
32.8 m in total position, i.e., smaller (about half) than the 
short arc overlap values, giving us confidence that the longer 
arcs are appropriate for building normal equations. 

Normal equations of the historical lunar radiometric track- 
ing data, recently analyzed for the creation of the GLGM-3 
gravity field (Mazarico et al. 2010), were recomputed with 
GLGM-3 as an a priori. They were then combined with the 
LRO normal equations and the gravity field inversion was 
performed. Even though our a priori weighting of the radio- 
metric data was appropriate, the altimetric crossover data had 
to be downweighted, because of the large-amplitude oscilla- 
tions they induced in the farside gravity field, even with the 
presence of a Kaula constraint. We also needed to discard 
the crossover normal equations of a few phases, before the 
obtained farside gravity field did not show excessive power 
(compared to GLGM-3 and SGMIOOi from Goossens et al. 
2011). While we continue to tackle this issue, a possible 


Table 3 Average overlap RMS differences, in meters, in the along-track (A), cross-track (C) and radial (R) directions, as well as in total position (T) 


A priori gravity 

Radio- 

-only orbit overlaps (m) 


Radio+crossover orbit overlaps (m) 



Along 

Cross 

Radial 

Total 

Along 

Cross 

Radial 

Total 

GLGM-3 all phases 

47.57 

49.91 

4.50 

70.06 

17.22 

14.26 

3.96 

22.91 

GLGM-3 NO_01 to NO_13 

41.81 

40.87 

3.47 

58.77 

16.38 

13.72 

3.66 

21.82 

LLGM- 1 all phases 

16.71 

15.44 

1.58 

23.39 

10.18 

8.37 

1.80 

13.63 

LLGM-1 NO_01 to NO_13 

10.67 

10.09 

1.13 

14.79 

8.35 

8.21 

1.62 

12.00 


For each a priori gravity field (GLGM-3 or LLGM-1), the averages are calculated over all 21 phases and over the 13 nominal phases. Results with 
and without the altimetric crossovers are given. More details (phase averages) are given in Tables S2 and S3 
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explanation is the sparsity of the crossover data of the farside 
(Fig. 2c-d): some data can outweigh the Kaula constraint, but 
are too few to yield robust estimates of the spherical harmon- 
ics expansion coefficients. 

Our selected gravity field, LLGM- 1 (LRO Lunar Gravity 
Model), is not manifestly superior than GLGM-3 in terms of 
geophysical performance, based on admittance and correla- 
tion studies (following Mazarico et al. 2010; not presented 
here). Coefficient differences between LLGM- 1 and GLGM- 
3 are shown in Figure S5. The low-order coefficients at high 
degree show stronger changes compared to other high-degree 
coefficients. This is consistent with the LRO data affect- 
ing the determination of low-order coefficients through the 
m-daily perturbations. The relatively large coefficient differ- 
ences in the / = 5-70 band are due to a combination of the 
Kaula-type decrease in coefficient power with degree, and to 
the re-arrangement of the lumped coefficients, affected by the 
farside gravity field. When expanding both the GLGM-3 and 
LLGM- 1 fields, we find that the gravity anomaly changes are 
small on the nearside (<10mgal RMS excluding the polar 
regions) and larger on the farside, as expected (~62 mgal 
RMS). We did expect a better determination of the farside 
gravity anomalies thanks to the altimetric crossovers. How- 
ever, in the non-polar regions, as mentioned above, the cov- 
erage by LOLA crossovers actually used in the solution is 
very sparse. In the polar regions (down to ~70° latitude), 
the lunar farside gravity field was already well character- 
ized by the historical tracking data. However, due the great 
majority of LOLA crossovers occurring in that region, we 
see larger changes than on the nearside (~42 mgal RMS). 
Those are appreciable, although significantly smaller than 
the ~83 mgal RMS observed in the deep farside (spherical 
cap of half-angle 60°, centered on 180°E, 0°N). Multi-month 
arcs can allow us to access the complete set and potentially 
improve the farside gravity anomalies. This computationally 
intensive approach may be the subject of future work. Our 
goal here was to obtain a gravity solution “tuned” to LRO, 
with the orbital performance as metric. 

Nevertheless, the addition of the LRO dataset can poten- 
tially benefit the estimation of low-degree terms in the gravi- 
tational signal, given the longer arc duration (5 days) and their 
links, through crossovers, over monthly timescales. Table 2 
lists our estimates in comparison to previous results (Dickey 
et al. 1994; Williams et al. 2009, 2011; Konopliv et al. 
200 1 ; Goosens and Matsumoto 2008 ; Matsumoto et al. 20 1 0; 
Goossens et al. 2011; Mazarico et al. 2010). The LLGM-1 
adjusted values are generally consistent with the existing 
solutions, in particular LP100K and LP150Q (Konopliv et al. 
2000, 2001), except for C 21 and kj. Our C '21 estimate is rather 
large and of negative sign, compared to previous estimates 
that were generally small and positive. Although it is not 
correlated with .S '2 1 (—0.07), the covariance matrix shows 
that the correlation coefficients of C 21 with C 30 , C 32 and 


C 41 are rather large (0.64, —0.79 and —0.65, respectively). 
They are, however, consistently smaller than in the GLGM-3 
solution, likely due to the addition of the LRO data. Our k .2 
value is very close to our a priori value based on Williams 
et al. (1987). Other groups’ estimates are consistently lower, 
in the 0.023-0.026 range, as is our GLGM-3 adjustment 
(0.0235 ± 0.0002). We note however the latest SELENE 
estimate using VLBI (Goossens et al. 2011) data is clos- 
est to our number (0.0255 ± 0.00016). Estimates based on 
Lunar Laser Ranging, after a long decrease due to model- 
ing change over the past decade (0.0302 by Dickey et al. 
1994 to 0.0199 by Williams et al. 2009), have recently been 
revised to values more consistent with satellite-based results 
(0.0229 ± 0.0020 by Williams et al. 2011). The correlation 
of our k .2 estimate with the other low-degree coefficients is 
small, with a maximum (absolute) value near 0.27 with .S' 30 , 
followed by 521 , C 22 and S 22 (smaller than 0.25). This is an 
improvement over GLGM-3, where the maximum correla- 
tion was 0.3, with 52 1 • 

3.4 Precise orbit reconstruction 

To further improve the orbit accuracy, we reprocessed the 
LRO tracking data, by using the new LLGM- 1 gravity field 
as an a priori to reprocess the short arcs. The radiometric- 
only arcs show improved residual Doppler and range RMS, 
with maxima around 1 mm s _ 1 and 4 m, respectively. This 
is mainly due to a better fit when the orbit is seen edge-on 
(Fig. 3), presumably because of a better (LRO-tuned) mod- 
eling of the integrated farside gravitational effects. 

The radiometric-only orbit performance is also excellent, 
with dramatic improvements compared to GLGM-3 in all 
directions (see Fig. 4 and Table S4). The mission-averaged 
total position knowledge is around ~23 m RMS (Table 3), at 
the level obtained with altimetric crossovers using GLGM-3 
as the a priori model. Moreover, we note that if we take the 
average performance over the nominal mission only, the orbit 
quality is further improved to ~ 15 m RMS. While this could 
be due to the specifics of the later phases (SM_01 to SM_05), 
it is plausible that those latter phases would see additional 
improvements with a gravity solution that includes data from 
those phases. In other words, the overlap numbers probably 
reflect orbit reconstruction performance for the earlier phases 
and something closer to orbit prediction performance for the 
latter ones. This will be answered by continuing work and 
iterations on the LRO gravity fields. 

With the inclusion of the altimetric data in the orbit 
convergence, similar observations can be made, but the orbit 
consistency is strengthened even further (Table S4). The 
calculated orbit overlap RMS over the whole mission are 
~8.3, 8.2 and 1 .6 m in the along-track, cross-track and radial 
directions, respectively, and 1 2.0 m in total position. The rela- 
tive improvement is much smaller than the 70% obtained with 
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Fig. 6 Topographic maps in polar stereographic projection of the lunar north pole region (axis are in kilometers), created from LOLA data 
geolocated with the navigation orbits (a) and with the LLGM-1 orbits (b). The map resolution is 20 m (true scale at the pole) 


GLGM-3, but still significant (19%) given the magnitude of 
the radiometric-only errors. The gain in the radial direction 
is larger (about 50%), but does not yield significantly smaller 
crossover discrepancies (Table 1). We note that on average 
the radial direction is slightly degraded by the addition of the 
altimetric crossovers. However, the crossover discrepancies 
are now commensurate with the radial overlap RMS values. 
The crossover residuals are smaller near the poles, where the 
majority of crossovers occur. The GLGM-3 overlaps were 
too large to observe this effect, but it was expected from the 
Rowlands et al. (2009) LRO simulation study: the crossovers 
tend to distribute the orbital errors among all directions, cor- 
rupting the empirical rule of the radial errors being ten times 
smaller than the horizontal errors. That rule still holds rela- 
tively well for the radiometric-only (both with GLGM-3 and 
LLGM-1), leading to larger radial errors when including the 
LOLA data. We find that the orbit differences between the 
radiometric-only orbits and those that include the crossovers 
(10.8, 7.3, 1.5 and 14.2 m on average in the ACR directions 
and total position, respectively) are commensurate with the 
differences in their overlap RMS values (Table 3), which 
is consistent with the crossovers bringing largely favorable 
orbital changes. 


4 Discussion and closing remarks 

4. 1 Orbit accuracy assessment from independent data 

The overlap analysis above shows that the LRO orbits are 
consistent at the 10- to 15-m level. This is technically a 


measure of precision and not of the orbit accuracy itself 
(as required for the LRO position knowledge; Vondrak et al. 
2010). The orbit accuracy is difficult to assess with the track- 
ing measurement types we used, because the true trajectory is 
of course unknown. Here, we use independent data to show 
that the orbit precision numbers (Sects. 3.2, 3.4) are com- 
mensurate with the orbit accuracy. 

Although we used a small amount of LOLA data in the 
POD process to construct the altimetric crossovers, the over- 
whelming majority of the altimetric measurements (laser 
bounce points) did not directly contribute to the orbit adjust- 
ment. The topographic maps produced from the complete 
LOLA dataset can thus be viewed as a legitimate test of the 
orbit accuracy, especially in the polar regions where the data 
density is high. In polar stereographic projection, with the 
orbits rotating with a monthly period, the sharpness of the 
topographic maps is possible only with accurate orbits. 

In Fig. 6, we compare maps of the lunar north polar region 
(89-90° N) at 20-m resolution created with the navigation 
orbits (Fig. 6a) and with the orbits presented in Sect. 3.4 
(Fig. 6b). Although artifacts due to orbit errors are still pres- 
ent in Fig. 6b, they are significantly reduced and are the result 
of a small number of tracks. Overall, the geolocation is excel- 
lent and broadly consistent with the overlap RMS decrease 
compared to other solutions (Table 3), and with the expected 
20-m precision level determined from overlaps. 

Images acquired with the LROC NAC camera, at 50-cm 
resolution, can also provide strong measures of the orbit 
accuracy. The LROC instrument is repeatedly pointed at the 
Apollo landing sites, as Constellation sites, but also as opti- 
mal calibration targets. To calibrate the instrument pointing 
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Table 4 Summary of the 
statistics of the timing residuals 
within LROC NAC images of 
the predicted and observed 
positions of Apollo/Lunokhod 
artifacts 


Results with both the navigation 
and the LLGM- 1 orbit 
reconstructions are given 


Target 

Number of 
observations 

Timing errors (ms) 
navigation orbits 

Timing errors (ms) 
LLGM-1 orbits 

Mean 

RMS 

Mean 

RMS 

Apollo 11 

10 

-7.19 

11.12 

-2.83 

6.19 

Apollo 12 

11 

-9.98 

28.58 

0.44 

6.13 

Apollo 14 

13 

-6.36 

12.31 

-0.54 

12.36 

Apollo 15 

11 

-4.75 

7.18 

-5.38 

6.56 

Apollo 16 

10 

-5.16 

19.61 

-4.62 

12.84 

Apollo 17 

11 

-18.46 

36.54 

-4.40 

8.96 

Lunokhod 1 

7 

-20.52 

35.21 

-2.98 

9.3 

Average 

10.4 

-10.35 

21.51 

-2.90 

8.91 

Weighted average 

- 

-9.79 

20.67 

-2.81 

8.96 


and timing, the LROC team has been keeping track of the 
position of known human artifacts in the LROC NAC image 
frames. Those data (M.S. Robinson, personal communica- 
tion) show that the orbits obtained in this work achieve 
improved accuracy, reducing the scatter of timing residu- 
als (between expected and actual) from ~21 ms RMS to 
9 ms RMS (Table 4, Figure S6). The mean of the residuals 
decreases significantly as well, from about —10 to —3 ms. 
Given the ~ 1.7 km s -1 LRO orbital velocity, those numbers 
are in good agreement with the results of the overlap analysis 
(i.e., ~15 m RMS, with an average of ~5 m). 

4.2 Creation of products and distribution 

The orbits presented in Sect. 3.4 are made available as part 
of the LRO Radio Science PDS archive (http://imbrium.mit. 
edu/LRORS/). They are distributed as a set of SPICE ker- 
nels (Acton 1996), each spanning one of the monthly phases 
described in Sect. 2.1. A linear taper is used in the overlap 
periods to combine the trajectories of consecutive arcs. The 
LLGM- 1 gravity field is also made available through the PDS 
archive or by contacting the authors. 

4.3 Summary of results 

In this paper, we presented the data and methods used by 
the LOLA POD team to reconstruct the LRO spacecraft tra- 
jectory to satisfy the LRO position knowledge requirements 
and to enable the geolocation and comparison of the various 
high-resolution datasets being acquired by the instruments 
onboard the spacecraft. We demonstrated through the use of 
overlap analysis and with the help of additional independent 
datasets that the orbits obtained with radiometric tracking 
data supplemented by LOLA-derived altimetric constraints 
achieve position accuracy levels around 20 m in total position 


and near 1 m radially. A gravity field solution that combined 
pre-LRO tracking data with the newly acquired LRO geodetic 
measurements was developed and tuned to further improve 
the LRO trajectory. This resulted in accuracies better than 
20 m over the whole mission (July 2009 to January 2011). 
This work will be continued to provide a consistent accuracy 
of the trajectory for all the data acquired by LRO and to estab- 
lish a high-quality lunar reference frame to combine datasets 
of multiple spacecraft, past or future. Further enhancements 
can also be expected: by the creation of multi-month arcs to 
make use of a larger number of altimetric crossovers; and ulti- 
mately by a reprocessing of the data with a much-improved 
gravity held solution obtained by the GRAIL mission (Zuber 
etal. 2011). 
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